Mon. Not. R. Astron. Soc. 000, HUTU (20081 Printed 3 September 2008 (MN M£X style file v2.2) 



OO 
O 



00 



Microstructure and kinematics of H2O masers in the 
massive star forming region IRAS 06061+2151 

K. Motogi, 1 * Y. Watanabe, 1 K. Sorai, 1 ' 2 A. Habe, 1 ' 2 M. Honma, 3 ' 4 H. Imai, 5 ' 6 
A. Yamauchi, 7 ' 8 H. Kobayashi, 3 ' 4 K. Fujisawa, 9 T. Omodaka, 5 ' 6 H. Takaba, 10 
O : K. M. Shibata, 3 ' 4 T. Minamidani, 1 ' 2 K. Wakamatsu, 10 H. Sudou, 10 
E. Kawai, 11 and Y. Koyama 11 

1 Department of Cosmosciences, Graduate School of Science, Hokkaido University, N10 W8, Sapporo 060-0810, Japan 
2 Department of Physics, Faculty of Science, Hokkaido University, N10 W8, Sapporo 060-0810, Japan 

3 Department of Astronomical Science, The Graduate University for Advanced Studies, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 
' 4 Mizusawa VERA Observatory, National Astronomical Observatory of Japan, 2-12 Hoshigaoka, Mizusawa, Oshu, Iwate 023-0861, Japan 
^Department of Physics, Faculty of Science, Kagoshima University, 1-21-35 Korimoto, Kagoshima 890-0065, Japan 
< ' i 6 Graduate School of Science and Engineering, Kagoshima University, 1-21-35 Korimoto, Kagoshima 890-0065, Japan 
i-^H ■ 7 Institute of Physics, University of Tsukuba, 1-1-1 Ten-nodai, Tsukuba, Ibaraki 305-8571, Japan 
Or 8 Mizusawa VERA Observatory, National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 
^ 9 Department of Physics, Faculty of Science, Yamaguchi University, 1677-1 Yoshida, Yamaguchi, Yamaguchi 753-8512, Japan 
» I , 10 Faculty of Engineering, Gifu University, 1-1, Yanagido, Gifu 501-1193, Japan 

' 11 Kashima Space Research Center, National Institute of Information and Communications Technology, 
^ 893-1 Hirai, Kashima, Ibaraki 314-8501, Japan 

Accepted 2008 July 14. Received 2008 July 14; in original form 2008 April 29 

in 

^ ; ABSTRACT 

■ We have made multi-epoch VLBI observations of H2O maser emission in the mas- 

sive star forming region IRAS 06061+2151 with the Japanese VLBI network (JVN) 
from 2005 May to 2007 October. The detected maser features are distributed within 
an l"xl" (2000 aux 2000 au at the source position) around the ultra-compact Hll 
I region seen in radio continuum emission. Their bipolar morphology and expanding 

motion traced through their relative proper motions indicate that they are excited 
by an energetic bipolar outflow. Our three-dimensional model fitting has shown that 
the maser kinematical structure in IRAS 06061+2151 is able to be explained by a 
biconical outflow with a large opening angle (> 50°). The position angle of the flow 
major axis coincides very well with that of the large scale jet seen in 2.1 /jm hydrogen 
emission. This maser geometry indicates the existence of dual structures composed of 
a collimated jet and a less collimated massive molecular flow. We have also detected a 
large velocity gradient in the southern maser group. This can be explained by a very 
small (on a scale of several tens of au) and clumpy (the density contrast by an order of 
magnitude or more) structure of the parental cloud. Such a structure may be formed 
by strong instability of shock front or splitting of high density core. 

Key words: ISM: jets and outflows - masers - stars: formation - ISM: kinematics 
and dynamics. 
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1 INTRODUCTION 

The studies of massive star formation are essentially impor- 
tant for astrophysics. Because of their powerful protostel- 
lar or stellar winds and radiation, their formation is deeply 
conce rned with activities of lo cal star formation in galaxies 
fe.g.. lElmegreen fc Ladalll977f ). Once massive young stellar 
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objects (MYSOs) are formed, their strong activities imme- 
diately interact with surrounding clouds. This interaction 
often results in fragmentation and dispersion of the parental 
cloud, and sometimes stops following star formation. How- 
ever, another star formation is triggered frequently in such 
compressed fragment. In spi te of numerous observa tions and 
numerical simulations (e.g.. lElmegreen et al. I ll998l . and ref- 
erences therein), quantitative and detailed studies on such 
interaction are still insufficient. 
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Direct observations of local kinematics around MYSOs 
would give intrinsic decision for such studies. If detailed 
information of the motion of gas in these interactions is 
available, we can examine the triggering effect quantita- 
tively. However, observational studies of their kinematics 
have some serious difficulties. First, all of MYSOs are still 
deeply embedded in the parental molecular clouds because 
of their rapid evolution. Second, they are usually seen in 
crowded star clusters. Third, they are far distant from the 
Sun (>500pc). In such deep and crowded massive star form- 
ing regions, star formation process is complex since several 
effects, such as stellar winds, radiations and stellar or inter- 
stellar magnetic fields, etc, are mixed. Consequently, high 
angular resolution is essential to resolve into individual ef- 
fects and local kinematics. 

Observations of H2O maser emission with a very long 
baseline interferometre (VLBI) is very suitable for such stud- 
ies. Because the maser excitation requires a high density 
(riH 2 ~ 10 9 cm -3 ) and a temperature (~600K), H2O maser 
emission in a massive star forming region is usually asso- 
ciated with strong shock front induced by MYSOs (e.g., 
lElitzudll992h . Each maser source is composed of small gas 
clumps called maser feature, which has typically a size of 
1 au and line of sigh t velocity width of a few km s _ (e.g., 
iReid fc Moranlll98ll ). This scale is enough to resolve the spa- 
tial and kinematical structure mentioned above. Moreover, 
extremely high spatial resolution of VLBI obsevations enable 
us to detect small proper motions of a few milli-arcseconds 
(mas) per year, which reflects the local kinematics pro- 
jected onto the celestial plane. For these reasons, H2O maser 
emission is a good tracer of three-dimensional kinemat- 
ics around MYSOs. In fact, many VLBI observations have 
shown highly systematic a nd directional proper motions 
of H2O maser features (e.g., Genzel et ah! 198ll; Imai et all 
200d; " 
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Moscadelli et al.ll2007l ). which are thought to trace propa- 



gation of shock fronts induced by a bipolar outflow, a pro- 
tosteller wind, or an expansion of ultra-compact (UC)HII 
region. 

In this paper, we report results of the VLBI observations 
of H2O maser source WB755 in the massive star forming 
region IRAS 06061+2151. IRAS 06061+2151 is one of the 
most mas sive star forming region in the Gemi ni OBI cloud 
complex (|Carpenter, Snell fc Schloerbl 1995bT). Its photo- 
metric distance is 2.0 kpc (jCarpenter. Snell fc Schloerbl 
Il995al ). Several observations h ave indicated active star 
form ati on in this region (e.g., Hanson. L uhman & Rickc 
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Thompson et al.ll2006l ). This region is deeply embedded in 
the extended (~1 pc) and massiv e (^7000 Mq) m olecular 
cloud complex observed in C ls O dSaito et~al 1 l2007h . There 
is a compact clus ter seen in near-infrare d (NIR), AFGL 
5182 at the centre ^Carpenter et al.lll995bl). Five bright Ks- 
band sources are identified in this cluster ( Anandarao et aD 
120041 ). They have c oncluded these sources are early B-type 
protostars (S1-S5 in l Anandarao et al.ll2004T ). They also have 
detected knot-like structure which suggests a protostellar 
jet. In addition, three weak centimetre continuum sources 
have been detected with the Very Large Array (VLA) in B 
configur ation, and all of them are though t to be UCHII 
regions ([Kurtz. Churchwell fc Woodlll994h . Only the most 
southern and brightest UCHII region is coincident with 



the sourthern Ks-b&nd source, S5. The other two UCHII 
regions are also very close to S4 source, but they are 
clearly individual sources. The geometry between these five 
_R"s-band sources and three cent i metre continuum sources 
is clearly described in lBik et al.l (|2006). The maser source, 
WB755, is associated wit h the most north ern UCHII 
region, G188.794+1.031, in iKurtz et""aH (|l994T ) (see section 
3). This region is too small to resolve with 1" resolution in 
the VLA observation at 8.4 GHz. This scale corresponds to 
about 2000 au at 2.0 kpc and abou t ten times smalle r than 
the typical size of UCHII region l|Hoare et al.ll2007T ). It is 
probable that this UCHII region is extremely young and 
its ionization source is still in the accretion phase. 

We describe the observations and data reduction in sec- 
tion 2, and the observational result in section 3. Discussions 
and model fitting are described in section 4. We summarize 
our conclusions in section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

Observatio ns were made w ith the Japanese VLBI Network 
(JVN, e.g- lDoi et ai1l2006r i at six epochs from 2004 May to 
2007 October in H 2 (Jif a x b = 6i6-5 23 ) at 22.23508 GHz. 
The baselines of the JVN is ranged from 50 to 2560 km, and 
the maximum angular resolution is about 1 mas at 22 GHz. 
This resolution corresponds to 2.0 au at 2.0 kpc. 

WB755 was observed for 8-10 hrs in each epoch. Other 
two continuum sources, J0530+13 and OJ287, were also ob- 
served in every 30-60 minutes for phase and bandpass cali- 
brations. The received signals were recorded with the SONY 
DIR1000 recorder. Their sampling rate is 128 Mbps with 
2-bit sampling. The d ata correlation was made with the 
Mitaka FX correlator (|Shibata et al.lll998T l. The correlated 
data are divided into 1024 spectral channels. Each chan- 
nel have a velocity resolution of 0.21 kms^ 1 and full band- 
width is 16 MHz which corresponds to a velocity coverage 
of 215 kms -1 at the frequency. Table 1 summarizes each of 
the observations. 

Data reduction was made using the National Radio As- 
tronomy Observatory (NRAO) Astronomical Imaging Pro- 
cessing System (AIPS). Amplitude calibration of visibility 
was made with the template spectru m method with the to - 
tal power of strong maser emission (jDiamond et al.lll995T ) 
in the AIPS task ACFIT. Instrumental delays and band- 
pass characteristic were removed using the data of the two 
continuum calibrators in the task FRING and BPASS. Un- 
certainty in the estimated delays was within 5-7 ns. Be- 
cause the detected maser emission is found within 1.2 MHz 
from the phase- reference channel, this error corresponds to 
the error of relative position within 30 micro arcseconds 
(/xas) for a baseline of 1000 km. We have selected the ve- 
locity channel of intensity peak for phase-reference chan- 
nel in each epoch. A Doppler velocity of referenced channel 
is around -8.0 km s _1 in all the epochs. Fringe- fitting and 
self-calibration procedure were made for the reference chan- 
nel and their solutions were adopted for other channels of 
same epoch. Because the reference channel often has several 
maser components, the self-calibration was made by hybrid 
mapping method (e.g.. iReadhead fc Wilkinsonlll978h . After 
these calibrations, we have searched maser components in 
the 1.5" x 1.5" area of each channel map. Finally, we have 
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Table 1. Detailed status of each epochs 



Epoch 


Duration 


Antennas" 


Reference 6 


1-cr 


Synthesized 


PA 








Velocity (LSR) 


noise level 


Beam 






(hr) 




(kms" 1 ) 


(mjy beam -1 ) 


(mas X mas) 


(°) 


2005/5/13 


9 


MZ, IR, OG, IS 


-9.67 


90 


1.2x1.0 


-79 


2006/1/30 


9 


MZ, IR, OG, IS, KS, TM 


-7.77 


35 


1.5x1.0 


-71 


2006/4/17 


9 


MZ, IR, OG, IS, KS, TM 


-7.56 


23 


2.0x0.9 


-43 


2007/1/29 


8 


MZ, IR, IS, KS 


-8.14 


15 


2.5x0.8 


-46 


2007/5/1 


8 


MZ, IR, IS, OG, KS, TM 


-8.08 


63 


1.5x1.0 


-78 


2007/10/21 


10 


MZ, IR, OG, IS, KS 


-8.08 


38 


2.1x0.8 


-43 



a MZ, IR, OG, IS The VERA 20 m antennas at Mizusawa, Iriki, Ogasawara, Ishigaki. 

KS : The NICT 34 m antenna at Kashima. TM : The Hokkaido University 11 m antenna at Tomakomai. 
6 Phase-referenced velocity channel 



performed two-dimensional Gaussian fitting to estimate the 
position of brightness peak for detected maser components 
with the task SAD. Uncertainty of this fitting is typically 
30-70 ^as. 

The individual maser clump usually appears in suc- 
cessive velocity channels, since the maser emission has line 
width of a few km s _1 . The maser component in each ve- 
locity channel is commonly called a maser spot. We adopt a 
term of maser feature, which is a small gas clump as men- 
tioned above and appears in successive velocity channels, to 
form single-peaked emission in this paper. That is, a maser 
feature consists of several maser spots which are located 
within 2.0 au and 3.0 km s _1 from the brightness-peak spot 
in the case of IRAS 06061+2151. Brightness-peak channel 
in each feature is assumed to be as Doppler velocity of the 
maser feature, while the intensity- weighted average of the 
positions of the associating maser spots is adopted as its 
position. Detailed information of detected maser feature is 
given in Table 2. 



3 RESULTS 

3.1 Absolute coordinates 

First, we roughly estimated absolute position of phase- 
reference maser spot at the sixth epoch, which is the bright- 
est spot and belongs to S8 feature in table 2, through the 
fringe-rate analysis. The estimated position is (a,5)j2ooo.o 
= (06 h 09 m 07 s .0, +21°50'41".2), which is roughly coincident 
with the peak emission of the rad io continuum source found 
with the VLA l|Kurtz et al.ll 19941 ). Since the phase calibrator 
was scanned every 30-60 minutes, short-term fringe-phase 
fluctuation cannot be removed, and the estimated position 
has uncertainty of about 150mas. This error caused by resid- 
ual phase error is, however, enough small to be compared 
with the diametre of the continuum source (jKurtz et al.l 
1 19941 ). 

Because detected maser features are seemed to be com- 
pletely inside of unresolved radio emission (see next subsec- 
tion), a geometrical relationship between the maser emission 
and the UC H II region still has inevitable uncertainity, but 
this result indicates that WB755 must be excited by this un- 
resolved radio continuum source, although the maser source 
has been thought to be located at the m ore eastern position 
(i.e. IRAS position; [Sunada et ai"1l2007l ). 



3.2 Spatial and velocity distribution 

The distribution of all maser features is presented in Fig 1. 
Each point corresponds to an individual maser feature, and 
its colour represents the Doppler velocity with respect to the 
local standard of rest. The offset of each feature is measured 
from the position-reference maser feature at (a,<5)j20oo.o = 
(06 h 09 m 07 s .O, +21°50'41".8) for the first three epochs. An- 
other position-reference feature have been chosen for the 
later three epochs, since the former feature disappeared. 
These features were detected at the same time in the third 
epoch. They are crowded within lOau in space and 0.3kms~ 
in Doppler velocity, hence, their relative motion is thought 
to be also small during the short period. This enables us 
to assume that their relative position is conserved between 
third and forth epoch. 

Maser features have concentrated within an area of 
l"xl" at all the epochs with LSR velocities ranged from 
5.0 to -20.0 kms" 1 . They are divided into the northern and 
middle groups, whose velocities are red-shifted with respect 
to the systemic velocity of about -l.Okms -1 (e.g., Sne ll et all 
1988; Sai to et al.l2007h . and the southern blue-shifted group. 
This size of distribution is rather smaller than the upper 
limit to the size of the continuum source, which may suggest 
that these maser features exist on the boundary between the 
ionized and neutral regions. Generally, there is high density 
gas compressed by expansion of a ionized region in a such 
place, which should support maser excitation. 



3.3 Intensity and filamentary structure 

Intensity of each maser spot has ranged from several tenths 
to hundreds Jy beam -1 (Table 2). The intensity-peaked fea- 
ture is several hundred times brighter than the weakest one 
with the detection limit of 7 a. The intensity of each fea- 
ture also varies with epoch in a wide range (Table 3). Maser 
intensity mainly depends on the environment of excitation, 
the length of amplification and the direction of beaming, ft 
requires some geometrical model to explain these variations 
theoretically. 

Most of strong features, including the intensity-peak 
features mentioned above, are clustered in the southern 
group, and they exhibit filamentary alignment among all 
epochs. H2O maser emission is often excited in a dense gas 
layer behind shock front, since an ex tremely high density of 
10 9 cm" 3 and temperature of 600 K (|Elitzurlll992 ) are satis- 
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Table 2. Detected Maser Features 



Epoch 


Group" 


Vlsr 


Offset c ( 


mas ) 


Line Width 


Intensity 




& Number 6 


(km s ) 


X (err) 


Y (err) 


(km s ) 


(Jy beam -1 ) 


1 


Nl e 


4.35 


(0.02) 


(0.02) 


0.42 


1.77 




N2 


-0.71 


99.95 (0.02) 


-88.55 (0.02) 


0.03 


2.51 




C 1 

Si 


-5.7b 


i oi tt 1 o ni \ 

181.77 (0.01) 


a cr a to ( r\ r\ i \ 

-454. 7y (0.01) 


2.73 


20.1 




S2 


-5.55 


"i "7/" oo / o oo\ 

17b. 32 (0.02) 


A i~' O O "1 { C\ 00\ 

-4b8.81 (0.02) 


1.47 


0.32 




S3 


-5.55 


174.8 (0.0b) 


A TT OO { C\ O A \ 

-477.03 (0.04) 


0.03 


5.24 




S4 


-7.45 


1 A A C O 1 o m \ 

144. bo (0.01) 


trio oo I c\ r\ i \ 

-513.92 (0.01) 


1.08 


103 




S5 


-8.50 


173.02 (0.01) 


-502.92 (0.01) 


2.73 


149 




bb 


-9.77 


i cr t ot 1 o 01 \ 

157.37 (0.01) 


cr 1 o 11 / n r\t \ 

-513.11 (0.01) 


1.2b 


213 




S7 


-9.5b 


157.33 (0.01) 


r" 1 o t/^ /o ni \ 

-512.76 (0.01) 


1.05 


203 




b8 


-9.5b 


15/. 3 (0.01) 


trio qa / o 01 \ 
-512.54 (0.01 ) 


2.52 


204 




cn 

by 


-10. bl 


i cr i o rr I o c\-\\ 

151.05 (0.01) 


r 1 n cyf ( r\ Ol\ 
-513. OD (0.01) 


3.99 


14.8 




blO 


-14.40 


i cr o co 1 o o a\ 

150. b8 (0.04) 


r i n oo I o oo\ 

-512. 8y (0.03) 


0.21 


1.71 




C 1 1 

bit 


-19.4b 


i / \ rr /;■ i I r\ o 1 \ 

195. bl (0.01) 


A cr O TO I O OI \ 

-453.73 (0.01) 


1.08 


3.90 


2 


N3 


4.98 


15.75 (0.07) 


10.55 (0.05) 


0.42 


0.01 




Nl 


4.5b 


(0.12) 


(0.05) 


0.42 


0.90 




SI 


-5.55 


188.03 (0.07) 


-457.58 (0.02) 


1.48 


33.5 




S12 


-5.13 


189.98 (0.15) 


-474.72 (0.0b) 


0.21 


15.7 




S2 


-5.34 


188.77 (0.21) 


-474.b3 (0.04) 


0.42 


19.5 




S3 


-5.55 


188.49 (0.19) 


-483. iy (0.04) 


0.42 


18 




S13 


-5.55 


189.32 (0.11) 


-4bb.03 (0.03) 


0.42 


24.8 




S14 


-5.55 


1 OT o/" 1 1 o 1 A \ 

187.2b (0.14) 


ACT' OCT ( r\ oo\ 

-4bb.05 (0.08) 


0.21 


19.8 




c 1 n 

blo 


-5.55 


i oo r\ a l o i a \ 

192.04 (0.14) 


a c tr 1 o I o oo\ 

-4b5.18 (0.0y) 


0.21 


19.2 




bib 


-7.03 


i on ocr 1 o o cr \ 

13y.25 (0.05) 


rr o i oo I o o a \ 

-521.02 (0.04) 


0.42 


27.1 




bl7 


-7.bb 


inn c a 1 o o cr \ 

190. b4 (0.05) 


A<~\i o /o oo\ 

-493. b3 (0.03) 


1.00 


15b 




b5 


-8.08 


174.57 (0.02) 


ETO,i oo f o oi \ 

-504.32 (0.01) 


1.9 


424 




S8 


-8.5 


158.8 (0.03) 


n r O/l /o Ol\ 

-515.04 (0.01) 


2.32 


241 




Sb 


-8.71 


i cro cr /o no\ 

158.5b (0.02) 


trier oo f r\ oi\ 

-515.03 (0.01) 


1.9 


22b 




S18 


-9.5b 


137.8 (0.03) 


rr»i 1 o /o 

-521.18 (0.02) 


2.53 


122 




c 1 n 


-13.98 


i c\ oo /o oo\ 

lb0.92 (0.02) 


cror 1 ot ( o 01 \ 

-50b. 37 (0.01) 


3.17 


9.44 




S4 


-14.19 


"i a v if / o no\ 

145.15 (0.03) 


P"1T OO /O Ol\ 

-517.09 (0.01) 


4.22 


14.1 




S10 


-18.83 


"i T* o to /o i n\ 

150.73 (0.10) 


fir" r* /o nr\ 

-515.5 (0.05) 


0.21 


0.33 


3 


N4 


b.25 


5.87 (0.08) 


-5.78 (0.08) 


0.21 


0.18 






4.35 


5.48 (0.01) 


-5.09 (0.02) 


2.52 


0.78 




Nb 


4.14 


5.39 (0.05) 


-5.43 (0.05) 


0.21 


0.51 




Nl 


4.77 


(0.04) 


(0.05) 


0.03 


0.58 




Ml 


2.24 


238.91 (0.02) 


-295.71 (0.03) 


1.20 


1.5 




M2 


1.19 


245.98 (0.02) 


-30b.bl (0.02) 


1.05 


1.03 




M3 


1.19 


240.23 (0.05) 


-300.34 (0.05) 


0.21 


0.02 




S20 


-2.18 


133.73 (0.01) 


-442.52 (0.02) 


1.47 


22.2 




S21 


-4.29 


lib. 55 (0.05) 


-483.85 (0.04) 


0.42 


0.42 




SI 


-5.13 


188. bl (0.02) 


-457.89 (0.02) 


0.84 


0.29 




S22 


-4.92 


moo c i / f\ r\ a \ 

188.51 (0.04) 


/ir"o 10 /o r\ i ~>\ 

-458.18 (0.03) 


0.21 


4.91 




coo 

bid 


-7.24 


i a cr oo 1 o o a \ 

145. 2y (0.04) 


r 1 01 ( r\ a\ 

-518.21 (0.04) 


0.84 


34.4 




bo 


-7.87 


i Ttr no I o no\ 

175.02 (0.02) 


rr .1 Ti'ooiX 

-504.7 (0.01) 


2.52 


310 




S8 


-8.29 


"I P" o r" / o oo\ 

158.5b (0.02) 


r" 1 f to /o 

-515.72 (0.02) 


2.1 


125 




Sb 


-8.29 


"i c o p" o / o no\ 

158.53 (0.02) 


n c t /o ni \ 

-515.7 (0.01) 


1.89 


12b 




S7 


-9.13 


158.27 (0.0b) 


r - "i 00 /o 

-510.03 (0.05) 


0.21 


74. b 




CIO 

bl8 


-10.19 


i ot o i /o 01 \ 

137.91 (0.01) 


rr 1 to ( r\ oo\ 

-521.78 (0.02) 


1.20 


24.5 




cn 

by 


-11.03 


i cro ot 1 o no\ 

152.37 (0.03) 


cr 1 t i I c\ oo\ 

-517. ol (0.02) 


0.42 


10.4 




C/l 


- ±!.OD 


i/ic ik f'n rn'i 


no 10 /n c\\\ 


4.00 


10. y 




S10 


-17.35 


151.07 (0.01) 


-510.23 (0.02) 


1.47 


1.13 




S24 


-lb.93 


151.02 (0.04) 


-51b. 2b (0.04) 


0.21 


0.48 


4 


M4 


2.88 


240.01 (0.04) 


-295.38 (0.04) 


1.05 


0.85 




Ml 


2.45 


240.04 (0.03) 


-295.52 (0.04) 


1.20 


0.98 




M5 


2.24 


240.13 (0.0b) 


-295.y (0.0b) 


0.03 


1.12 




N5 


5.19 


5.48 (0.01) 


-5.09 (0.01) 


3.57 


11.2 




M2 


l.bl 


247.43 (0.02) 


-307.84 (0.02) 


2.1 


4.54 




N7 


2.45 


27.47 (0.11) 


20.9 (0.11) 


0.21 


0.31 




S20 


-2.b 


133. b2 (0.02) 


-443.02 (0.02) 


2.1 


23.1 




N8 


-2.81 


-18.52 (0.02) 


-97.bb (0.02) 


1.20 


14.0 




S25 


-5.7b 


125.35 (0.02) 


-529.38 (0.02) 


0.84 


7.92 
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Table 2 - continued 



Epoch Group" 


Vlsr 


Offset c 1 


mas ) 


Line Width 


Intensity 


& Number 6 


(km s ) 


X (err) 


Y (err) 


(km s ) 


(Jy beam -1 ) 


S26 


-5.76 


185.31 (0.09) 


-500.23 (0.09) 


0.42 


2.44 


S16 


-7.03 


138.87 (0.03) 


-524.46 (0.03) 


1.26 


15.2 


S27 


-6.82 


138.93 (0.02) 


-524.52 (0.02) 


0.84 


12.6 


S28 


-6.82 


137.41 (0.04) 


-526.88 (0.03) 


0.84 


6.57 


S5 


-7.45 


176.49 (0.03) 


-506.61 (0.03) 


1.47 


155 


S29 


-7.45 


176.64 (0.03) 


-506.53 (0.03) 


1.05 


153 


S30 


-8.29 


175.06 (0.07) 


-515.08 (0.07) 


0.42 


115 


S8 


-8.08 


158.56 (0.01) 


-518.32 (0.01) 


2.73 


275 


S31 


-10.61 


140.77 (0.02) 


-524.06 (0.02) 


1.68 


2.2 


OOO 

b23 


-10.61 


i a c oo 1 o oo\ 

145.83 (0.02) 


coi 01 t r\ oo\ 

-521.81 (U.U2j 


1.47 


3.19 


OOO 

b32 


-10.61 


i a o co 1 o oo\ 

142.53 (0.03) 


COI O A 1 c\ r\ A \ 

-521.84 (U.U4) 


1.05 


1.94 


ooo 

b33 


-10.61 


i a A o 1 o n c \ 

144.2 (0.05) 


C O A O/" 1 { C\ OC\ 

-524.26 (0.05) 


0.42 


1.99 


S18 


-11.45 


i o A ta 1 o oo\ 

134.74 (0.02) 


coc ot I c\ no\ 

-525.37 (0.02) 


1.26 


2.94 


S34 


-11.45 


i o o oo 1 o 1 n\ 

133.99 (0.10) 


coo o/^ /o r\n\ 

-533.86 (0.06) 


0.21 


1.53 


oo rr 

b35 


-12.3 


i co c a 1 o c\c\\ 

152.64 (0.09) 


-512.07 (0.06) 


0.21 


2.4 


S9 


-12.3 


"i c "i T o / o oo\ 

151.73 (0.03) 


coo co /o r\ i ~\\ 

-520.52 (0.03) 


1.05 


4.5 


S4 


-13.35 


i a n o 1 o no\ 

146.3 (0.02) 


coi oo I c\ no\ 

-521. 62 (U.U2) 


1.47 


2.32 


S36 


-15.67 


146.17 (0.02) 


CO"l O "1 /o 

-521.31 (0.02) 


1.68 


2.82 


b3/ 


-15.03 


i a a. nc; in nn\ 
14b. Uo (U.Uy) 


-521. ob (o.oy ) 


n oi 
U.21 


0.66 


S10 


-16.09 


i co T 1 o oo\ 

152.7 (0.03) 


c "i o or ( r\ oo\ 

-519.86 (0.02) 


0.84 


2.14 


5 N9 


6.67 


23.01 (0.12) 


6.56 (0.06) 


0.21 


0.95 


N10 


5.4 


23.08 (0.13) 


6.77 (0.06) 


0.21 


2.12 


N5 


4.77 


5.48 (0.04) 


-5.09 (0.02) 


2.94 


12.2 


M2 


0.98 


247.76 (0.15) 


-308.49 (0.08) 


0.21 


1.29 


M6 


0.77 


265.29 (0.12) 


-296.53 (0.06) 


0.21 


1.31 


S38 


-5.76 


ooo a o 1 r\ ir\\ 

203.48 (0.30) 


-488.47 (0.14) 


0.21 


4.01 


bo 


-8.08 


176.74 (0.10) 


COT OO I O O0\ 

-507.02 (0.03) 


1.05 


181 


S8 


-8.71 


158.84 (0.07) 


C"io o a i o r\ i ~>\ 

-519.04 (0.03) 


1.47 


53.6 


S35 


-10.4 


i co / ■ n l r\ c\c\\ 

152.69 (0.09) 


ci a r\c i c\ oc\ 

-514.06 (0.05) 


0.21 


1.75 


S39 


-14.4 


164.29 (0.06) 


C "1 o oo { c\ oo\ 

-510.02 (0.02) 


1.89 


3.12 


S3 


-14.19 


"1 OO / ■ i /o oo\ 

192.61 (0.39) 


A OO O A { O "1 A \ 

-490.84 (0.14) 


0.21 


1.15 


S36 


-14.4 


"1 AC* *7 /O O C \ 

146.7 (0.05) 


coo /o oo\ 

-522 (0.03) 


1.05 


3.15 


S37 


-14.4 


"i An i~' o /o c\r'\ 

146.68 (0.06) 


coo oo /o r\ i ~\\ 

-522.03 (0.03) 


0.84 


3.13 


S2 


-15.03 


~i r\ a 0"i /o /io\ 

194.01 (0.40) 


/lOO oo /o io\ 

-482.38 (0.13) 


0.21 


0.69 


S10 


-15.67 


"ICO CO /O "1 r?\ 

152.52 (0.17) 


coo t /o r\T\ 

-520.7 (0.07) 


0.21 


1.25 


S40 


-15.67 


"1 S ' O OO /O "1 o\ 

169.99 (0.13) 


COO TO { C\ O A \ 

-508.73 (0.04) 


0.84 


1.29 


6 N5 


4.56 


5.48 (0.01) 


-5.09 (0.01) 


2.52 


7.59 


N6 


4.35 


5.63 (0.03) 


-5.34 (0.03) 


0.42 


4.59 


M2 


1.61 


248.93 (0.02) 


-308.59 (0.02) 


2.1 


5.67 


M3 


1.82 


249.01 (0.03) 


-308.52 (0.04) 


0.63 


4.28 


S41 


-1.13 


121.22 (0.02) 


-530.02 (0.02) 


1.47 


1.01 


S42 


-4.08 


125.18 (0.02) 


-531.04 (0.02) 


1.26 


1.83 


S21 


-4.71 


115.82 (0.03) 


-484.89 (0.03) 


0.42 


0.9 


N8 


-3.03 


-17.89 (0.03) 


-98.08 (0.03) 


0.63 


1.27 


Nil 


-4.5 


35.03 (0.04) 


-148.2 (0.04) 


0.42 


0.75 


N12 


-4.5 


36.45 (0.03) 


-143.63 (0.04) 


0.42 


0.81 


N13 


-4.71 


35.1 (0.03) 


-148.37 (0.04) 


0.42 


0.98 


S43 


-5.76 


186.96 (0.02) 


-501.94 (0.02) 


0.84 


4.79 


S26 


-5.97 


187 (0.01) 


-501.78 (0.02) 


1.68 


6.67 


S44 


-5.97 


186.98 (0.02) 


-501.64 (0.02) 


1.26 


6.6 


S25 


-6.19 


122.83 (0.03) 


-531.84 (0.03) 


1.47 


3.63 


S23 


-7.45 


147.51 (0.03) 


-523.59 (0.03) 


0.84 


72.5 


S5 


-7.45 


178.33 (0.03) 


-508.45 (0.04) 


0.84 


40.3 


S29 


-7.24 


178.4 (0.03) 


-508.36 (0.03) 


0.63 


30.4 


S45 


-7.03 


178.4 (0.04) 


-508.38 (0.04) 


0.21 


13.4 


S8 


-8.08 


160.06 (0.01) 


-520.76 (0.01) 


1.89 


577 


S7 


-8.29 


159.9 (0.02) 


-521.04 (0.02) 


1.26 


446 


S33 


-9.77 


145.68 (0.02) 


-525.98 (0.02) 


0.63 


6.18 


S46 


-9.56 


145.78 (0.04) 


-525.68 (0.04) 


0.21 


4.93 


S47 


-9.56 


145.78 (0.03) 


-526.04 (0.04) 


0.21 


7.17 


S48 


-9.56 


144.79 (0.19) 


-534.37 (0.10) 


0.21 


2.38 


S4 


-9.98 


148.25 (0.21) 


-525 (0.32) 


0.21 


1.35 
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Table 2 - continued 



Epoch 


Group" 


Vlsr 


Offset c ( 


mas d ) 


Line Width 


Intensity 




& Number 6 


(km s -1 ) 


X (err) 


Y (err) 


(km s _1 ) 


(Jy beam ) 




S10 


-10.19 


150.96 (0.02) 


-524.46 (0.02) 


1.05 


2.48 




S18 


-11.66 


136.67 (0.01) 


-529.5 (0.01) 


1.89 


4.61 




S36 


-14.61 


147.59 (0.01) 


-524.11 (0.01) 


2.73 


4.8 




S37 


-14.4 


147.73 (0.02) 


-523.75 (0.02) 


0.63 


4.08 



a N: northern, M: middle, S: southern 

b Each feature, if it is detected, has same number among all the epochs. 

c Offset from the position-referenced feature. 

d Milli-arcsecond. 

e First position-referenced feature. 

f Second position-referenced feature. 
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200 


-200 




X [mas] 
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Figure 1. Destribution of all maser features detected in the 
six epochs. The coordinates indicates relative offsets from the 
position-referenced maser feature. Each point denotes each maser 
feature, and the colour indicates Doppler-velocity with respect 
to LSR. The location of the VLA 8.4 GHz continuum peak 
(G188. 794+1. 031) is indicated by the filled square, which has the 
uncertainty of 150mas because of the positional error of the phase- 
reference maser spot (see section 3.1). The location of MYSO 
estimated from our three-dimentional model fitting (see section 
4.2) is indicated by the filled triangle. The dotted and dot-dashed 
lines indicate the axis and opening angle of the best fitted bicon- 
ical outflow, respectively. 



fied for its excitation. Maser amplification also requires high 
velocity-coherence along its path. Consequently, maser am- 
plification occurs parallel to gas layer, and perpendicular to 
the shock propagation. The filamentary alignment of maser 
features, hence, indicates that there is a strong shock front 
around the southern region. 

Fig 2 shows a close-up view of the southern group. The 
features clearly exhibit arc-like alignment at the first and 
second epochs, which gradually changed to a straight-line 
with epoch. This indicates that there is a large gradient of 
the proper-motion velocity in the shocked gas layer. In this 
point of view, they seemed to consist of two different groups. 
The first group is the vertical part of the arc in the first 
epoch (V-wing) and the second one is the horizontal part 
(H-wing). It is noteworthy to mention that the intensity- 
peak maser is always located around the position where the 
two wings connect (i.e. the kink of the arc). This may also 
mean that the kink is most strongly shocked and results in 
extreme excitation of masers. 

3.4 Relative proper motions 

We calculate relative proper motions of maser features with 
respect to the position-referenced feature. Most of relative 
proper motions are larger than lmasyr -1 , corresponding to 
about 9.3 km s _1 at 2.0 kpc. These position drifts of maser 
features are much larger than the estimated error of relative 
position of ~0.1 mas. Details of the proper motions is given 
in Table 3. 

The relative proper motion vectors averaged for each 
feature are shown in Fig 3. All maser features exhibit ex- 
panding motion from the position-referenced feature (i.e. 
from the northern group). It is important that the relative 
proper motion vectors in the southern group are almost per- 
pendicular to their alignment. With mentioned above, maser 
alignment arises parallel to shock fron t and its motions of- 
ten reflect the shock propaga tion (e.g, iTorrelles et al.ll200ll ; 
llmai. Deguchi fc Sasa 0I I2OO2I ). If we assume that the maser 
features simply move along the shock propergation, these 
relative vectors are thought to be nearly same direction as 
the absolute proper motion vectors (i.e., NW-SE direction 
in this case.) 

If we also assume the symmetrical motion along this 
direction (i.e. like a bipolar flow), a half of relative motion 
should be good measure of absolute expanding motion. In 
this point of view, the largest relative proper motion faster 
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Table 3. All proper motions 



Group & Number 


Wing" 


Epochs 6 


Proper Motion c (mas yr 1 ) 
VX (err) VY (err) 


Velocity^ 
(km s _1 ) 


dV e LSR 
(km s _1 ) 


(Jy beam -1 ) 


N6 




3-6 


0.16 (0.04) 


0.06 (0.04) 


1.21 


0.21 


4.08 


N8 




4-6 


0.86 (0.05) 


-0.58 (0.05) 


7.38 


-0.21 


-13.3 


Ml 




3-4 


1.44 (0.05) 


0.24 (0.06) 


10.38 


0.21 


-0.517 


M2 




3- 4 

4- 5 

5- 6 


1.85 (0.04) 

l.O ^U.Ol ) 

2.46 (0.32) 


-1.56 (0.04) 
-0.21 (0.17) 


17.21 
20.92 
17.55 


0.42 
-0.63 
0.63 


3.51 
-3.25 
4.38 


Ml 
IVIO 




o o 


1.84 (0.04) 


-1.44 (0.04) 


1U.U-L 


u.uo 


o.uu 


SI 


V 


1- 2 

2- 3 


8.72 (0.10) 
2.77 (0.33) 


-3.89 (0.03) 
-1.47 (0.13) 


67.89 
22.30 


0.21 
0.42 


13.4 
-27.2 


S2 


V 


1- 2 

2- 5 


17.35 (0.30) 
a o f n ifi\ 

Q.Z I U.OO 1 


-8.12 (0.06) 

« 99 ff\ 1 1 *\ 


136.20 
53.36 


0.21 
-9.69 


13.2 
-18.8 


S3 


V 


1-2 
9 ^ 


19.08 (0.27) 
3.3 (0.35) 


-8.58 (0.09) 
-6.14 (0.12) 


148.74 




8 RA 


12.8 

1 Q 


S4 


H 


1-2 

3- 4 

4- 6 


0.72 (0.05) 
1.47 (0.03) 

9 (C\ 9CA 


-4.41 (0.02) 
-4 (0.03) 

K 07 (C\ AA~\ 


31.77 
30.30 
40.77 


-6.74 
-1.69 
3.37 


-88.9 
-11.6 
-0.97 


S5 


K 


1-2 

3- 4 

4- 5 

5- 6 


2.16 (0.03) 
1.87 (0.04) 
1 (0.41) 

Q Q/l (f\ 90\ 


-1.95 (0.02) 
-2.43 (0.04) 
-1.64 (0.17) 

-Q \ U. -L i- ) 


20.69 
21.80 
13.66 
31.92 


0.42 
0.42 
-0.63 
0.63 


275 
-155 
26 
-141 


S6 


H,K 


1-2 


l.DD ^U.U^t J 


9 fi£ (C\ ^9^ 


22.41 


1.05 


13 


S7 


H,K 


1-3 


1 (0.07) 


-3.52 (0.06) 


26.02 


0.42 


-188 


S8 


H,K 


1-2 

3- 4 

4- 5 

5- 6 


2.09 (0.04) 

U ^U.UO ) 

1.13 (0.29) 
2.569 (0.15) 


-3.06 (0.02) 
-o.o ^u.uo ^ 
-2.91 (0.13) 
-3.596 (0.07) 


26.35 
23.46 
22.20 
31.42 


1.05 
0.21 
-0.63 
0.632 


-23 
150 
-221 
523.4 


S9 


H 


1-3 


1.42 (0.03) 


-4.04 (0.03) 


30.45 


-0.42 


-4.4 


S10 


H 


1-2 

3- 4 

4- 5 

5- 6 


0.07 (0.15) 
2.07 (0.04) 
-0.75 (0.69) 
-3.26 (0.36) 


-3.63 (0.08) 
-4.62 (0.04) 
-3.37 (0.30) 
-7.89 (0.15) 


25.81 
35.99 
24.55 
60.70 


-4.42 
1.26 
0.42 
5.48 


-1.38 
1.01 
-0.89 
1.23 


S16 


H 


2-4 


-0.39 (0.06) 


-3.46 (0.05) 


24.76 





-11.9 


S18 


H 


3- 4 

4- 6 


-4.03 (0.03) 


-4.55 (0.03) 
-o.Dy ^u.uo ) 


43.22 
44.69 


-1.26 
-0.21 


-21.6 
1.67 


S19 


H 


3-4 


-U.ol (U.Uo J 


-o. / ) 


26.93 


-1.26 


-5.9 


S20 




3-4 


-0.13 (0.03) 


-0.64 (0.04) 


4.64 


-0.42 


0.9 


S21 




3-6 


-0.48 (0.04) 


-0.69 (0.03) 


5.98 


-0.42 


0.486 


S23 


H 


3- 4 

4- 6 


0.68 (0.05) 
2.31 (0.05) 


-4.58 (0.05) 
-2.46 (0.05) 


32.92 
23.99 


-3.37 
3.16 


-31.2 
69.3 


S25 


H 


4-6 


-3.47 (0.05) 


-3.4 (0.05) 


34.54 


-0.42 


-4.29 


S26 


V 


4-6 


2.32 (0.12) 


-2.13 (0.12) 


22.39 


-0.21 


4.23 


S29 


H 


4-6 


2.42 (0.06) 


-2.52 (0.06) 


24.84 


0.21 


-123 


S33 


H 


4-6 


2.03 (0.07) 


-2.36 (0.07) 


22.13 


0.84 


4.19 


S35 


H 


4-5 


0.17 (0.51) 


-7.98 (0.31) 


56.75 


1.9 


-0.65 


S36 


H 


4- 5 

5- 6 


2.11 (0.23) 
1.88 (0.11) 


-2.74 (0.14) 
-4.43 (0.06) 


24.59 
34.22 


1.26 
-0.21 


0.33 
1.65 
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Table 3 - continued 



Group & Number Wing" 


Epochs 6 


Proper Motion c (mas yr 1 ) 


Velocity d 


dV e LSR 








VX (err) VY (err) 


(km s -1 ) 


(km s _1 ) 


(Jy beam -1 ) 


S37 H 


4-5 


2.5 (0.44) -1.85 (0.38) 


22.11 


0.63 


2.47 




5-6 


2.2(0.14) -3.61(0.08) 


30.06 





0.95 



a H: H-wing, V: V-wing, K: Kink of arc. 

b Each two epochs indicates the duration in which proper motion is measured. 

c lmasyr - : 9.3 kms -1 at 2.0 kpc. 

d Absolute value of relative proper motion velocity. 

e Shift of Line of sight velocity in each feature between two epochs. 

f Shift of intensity in each feature between two epochs. 
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Figure 2. Close-up view of the southern group. Colour of the 
maser features indicates a detection epoch, epochl (blue), 2 (light 
blue), 3 (green), 4 (light green), 5 (yellow) and 6 (red). The kink 
area containing the brightest features is indicated by dashed cir- 
cle. 
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than 180 kms -1 is interpreted to the absolute proper motion 
of 90 km s . This value is rather fa ster than a typica l ex- 
panding velocity of UCHII regions l|Hoare et alj 120071 ). al- 
though such a bload line width of H2O maser e mission is seen 
in sev eral star forming regions (e.g., W49N; [g winnetaJj 
Il994al ) In addition to this, the velocity-coverage is about 
six times larger than the coverage of Doppler velocity. This 
difference i s possible in the c ase of a highly inclined out- 
flow (e.g., iGoddi et al. I 120051 ) rather than a spherically- 
symmetrical expanding flow. With these reasons, this ex- 
panding motion seems to trace a bipolar outflow rather than 
a spherical expansion of UC H II region. This is still indeter- 
minate at this point and our model fitting gives more deci- 
sion (See section 4.2). The dynamical time scale of the maser 
expansion estimated from the relative motions and the spa- 
tial distribution of about 1" is ranged from 10 2 to 10 3 yr. If 



Figure 3. Relative proper-motion vectors. An arrow indicates the 
feature-averaged relative proper-motion vector. The southern and 
middle groups show expanding motion from the northern group. 
In the southern group, a lot of relative motions are perpendicular 
to the maser alignment. 

their expanding motions trace a expansion of the ionization 
front, this small ionized region may be extremely young. 

A velocity gradient of the southern group is especially 
noteworthy. The V-wing shows remarkably fast motions (90- 
180 kms' 1 ) at the early epochs, while the H-wing shows 
relatively slow motions (30-65 kms -1 ). Of course, this large 
gradient clearly results in the change of their alignment men- 
tioned above. 

There exists steep inverse correlation between relative 
proper motion and intensity of the features (Fig 4). This 
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Figure 4. Relation between relative proper motion and inten- 
sity. Almost all of strong (> 30 Jy beam -1 ) features show small 
(< 5 masyr -1 ) motions. This steep inverse correlation can be 
interpreted as the effect of internal motion of gas layer. The best- 
fitting power law index is -1.6 and the correlation coefficient is 
0.42. 



is reasonable, since a feature wi th larger proper mot ion has 
larger internal velocity gradient (|Gwinn et al.lll994bl ), which 
prevents itself from obtaining coherent length. As a result, 
it may has shorter lifetime and weaker intensity. In fact, 
averaged intensity of the V-wing is relatively weaker than H- 
wing at all epochs, and the H-wing seems to be more stable 
than the V-wing which disappeared in the later epoch. 

Both wings show correlation between relative proper 
motion and position offset from the kink of the arc, and 
the V-wing has cleary steeper gradient than the H-wing in 
the proper motion-position offset diagram (Fig 5). In fact, 
several features of the V-wing show strong deceleration as 
they come down (see Table 3). This difference cannot be 
explained by an internal motion of single gas layer or a ge- 
ometry of typical gas flow swept-up with a simple bipolar 
outflow. 



4 DISCUSSION 

4.1 Interaction model 

These two wings show clearly different properties in their 
proper motion velocities and velocity gradient along the 
alignment. But their continuous alignment strongly suggests 
that they are excited in a single gas layer. With this point of 
view, we try to explain these differences with highly clumpy 
structure of their parental cloud. 



4-1.1 Motion of gas layer 

As a first step, we examine analytical treatment of motion 
of gas layer with simple plane parallel approximation, where 
a bipolar outflow strikes a uniform plane of parental cloud 
and sweeps up the gas into thin layer. In this compressed 
layer, one side is pressed by ram pressure and another side 
is supported by thermal pressure. An equation of motion for 
the gas layer is given by 




10 20 30 
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(b) V-wing 
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Figure 5. Velocity gradient in the both wings. The horizontal 
axis shows the one-dimentional (holizontal for H-wing and vertical 
for V-wing) position offset from the kink of arc (i.e. averaged 
position of S5 in table 2), and the vertical axis shows averaged 
proper motion in each position offset. Each solid line shows result 
of linear fitting (y = 0.04a; + 3.18 in the top; y = 0.20x + 4.29 
in the bottom). V-wing shows far steeper inverse correlation than 
H-wing. 



d ■ 2 

-j^(RpcloudR) = Pflowfflow! 



(1) 



where, p, v and R are the density, velocity and compressed 
length of the parental cloud, respectively. The left side cor- 
responds to time derivation of momentum per unit area of 
the layer. The right side of the equation corresponds to the 
ram pressure of the flow. 

We simply assume p c i ou d as constant. Then, equation 
(1) is transformed into 



R 



, "Jew 2 

v "flow 

^cloud 



RRY- 



(2) 



where, rifl ow and n c i ou( j are the number densities of molecular 
hydrogen in bipolar outflow and the parental cloud, respec- 
tively. R is still small in a short time after the interaction, 
hence we can assume R = and the equation (2) directly 
gives an initial velocity of the layer, 
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This means if the parental cloud has a density gradient, then 
single layer also has a gradient in an initial velocity. Such a 
velocity gradient simply affects maser intensity, since maser 
coherent path naturally increases in the slow part of gas 
layer, where gas density is likely to be dense, and maser exci- 
tation itself directly increases. Consequently, maser features 
in the slower part of gas layer is thought to show smaller 
proper-motions, stronger intensity and longer lifetime than 
the fast part. 

Because flow is not always perpendicular to gas layer, 
oblique shock wave is possible to occur. Such shock can bend 
the motion vector of gas layer, and these interactions can be 
directly detected as systematic change of radial velocity of 
features. In the case of IRAS 06061+2151, the V-wing seems 
to be more red-sifted than the H-wing at the early epochs, 
on the other hand, several high velocity features were decel- 
erated at the later epochs as mentioned above. Their line of 
sight velocities are blue-shifted comparable to the H-wing 
after the deceleration. This change of velocity vector can 
be interpreted as the effect of oblique shock. However, such 
deceleration itself is also possible to be explained by the in- 
crease of compressed length R, and the change of line of sight 
velocity might be an effect of internal motion of feature. In 
order to distinguish these effects, we need more sensitive 
observation enough to monitor faint part of thin gas layer. 

Molecular flow or protosteller wind is not a steady flow 
in practice. Some flows or winds have fluctuation, which 
causes time dependency of its ram pressure. It is well known 
that H2O maser emission in star for ming regions gener- 
ally has strong time variability (e.g., Mig enes et al.l Il999l ; 
iTrinidad et al.l 20031 ). It is possible that such variability is 
influenced by time dependence of ram pressure. The velocity 
of gas layer also depends on time in this case. If relaxation 
time of interaction is enough shorter than a cycle of fluctua- 
tion, a change of ram pressure immediately reflects in motion 
of compressed gas. We might be able to detect some corre- 
lation between systematic change of velocity and intensity 
in such a case. 

Fig 6 shows the schematic view of our interaction model. 
Denser gas clump is located at the southern side of cloud in 
this model. The kink of the arc corresponds to the densest 
part of clump, while density becomes lower with distance 
from the kink in the H-wing, as shown by proper-motion ve- 
locity gradient. On the other hand, the V-wing corresponds 
to diffuse part of cloud and faster than the H-wing, hence, 
the arc-like alignment is changed to straight-line. In this way, 
our model has given a reasonable explanation for change 
of alignment and proper-motion velocity gradient in the H- 
wing. If this qualitative explanation is correct, these behav- 
iors of maser gas layer tell us a distribution of a very small 
and dense gas clump surrounding a massive star forming re- 
gion. Since the difference of proper-motion velocity is up to 
5 times in the case of IRAS 06061+2151, corresponding den- 
sity gradient is estimated to be up to 25 times from equation 
(2). We discuss about this very clumpy structure in scale of 
several tens of au in the following subsection 4.1.2. 



4-1-2 Small density structure 

Maser excitation requires volume density of 10 9 cm -3 . This 
density is thought to be achieved by shock p ropagation. I n 
this case, pre-shock density is over 10 cm" 
which is slightly higher than that of typical molecular cloud 
core in scale of 0.1 pc. Such an extremely dense "bed" of 
maser is probably not intrinsic but induced by external ef- 
fect such as expansion of UCHII region. Even more, our re- 
sult means this "bed" is not uniform but extremely clumpy 
in scale of several tens of au. The lower limit of pre-shock 
density and the density ratio, which estimated from equa- 
tion (2), indicate that the density of clump is order of 10 7 - 
10 s cm -3 . 

There are two possibilities about the origin of such 
a dense and clumpy structure. One case is due to 
some external triggering. Previous simulations have shown 
that several types of instability of ionization front have 
been able to produce such dense and clumpy "bed" 
(e.g.. | Garcia- Segura fc Franco! Il99fj ; IWilliams et al. I Il999l ; 
iFranco et alj 120071 ). These instabilities arise as a result of 
expansion of UCHII region, and induce strong density gra- 
dient in surrounding neutral gas layer, which is thought to 
give effects on evolution of morp hology from UCH II re- 
gion to extended UCHII region (|Franco et al.ll2007f >. The 
most dense part becomes over ten times denser than the 
most diffuse part in these instabilities. Besides this, grav- 
itati onal instability also a rises in shocked gas layer it- 
self (jYoshida fc Habd I1992T ). All of these instabilities are 
thought to be mainly caused by anisot ropic pressure (ram , 
magnetic, etc) on a shock fro nt (e.g.. IWardle et all Il99fj| ; 
iGarci'a-Segura fc Francol 1 19961 ). However, typical scale of 
these simulations, which generally intend to study pc-scale 
structure, is not suitable for our observational scale smaller 
than 100 au. 

The other case is due to internal cause of self-gravity 
in the parental cloud core. If several high density parts are 
produced in the self-contracting cloud core at the same time 
on the formation of massive cluster, the parental cloud core 
have been intrinsically clumpy. Such high density clumps 
may be formed by splitting of high density core. Whether 
or not, there is no observational result which negates such a 
dense and small clump. 



4.2 Three-dimensional model fitting 

We fit the maser motions with a our three-dimensional 
model. From a bipolar-like distribution of maser fea- 
tures, we assume typ ical biconical outflow model (e.g., 
iMoscadelli et al] |2000| ). In this model, all maser features 
move along the surface of a cone. In addition, we also as- 
sume that the MYSO is located within the distribution of 
maser features. The free parametres in our model fitting are 
inclination of a cone from the line of sight, position angle on 
the celestial plane, opening ang le, and position of vertex (i.e. 
location of MYSO) as following iMoscadelli et ail l|2000l) . Be- 
cause of large uncertainty of position error among epochs, 
we also fit position of vertex as the free parametre for each 
epoch. From equation (1), the velocity of the gas layer that 
includes maser features is clearly related with the swept- 
up length R. Since the density of molecular cloud core in 
hydrostatic equilibrium generally shows power law distribu- 
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tion, we assume a velocity field on the surface of the cone 
under the constant or simple ram pressure, which decreases 
with R, as 

Mcone OC r", (4) 

where r is distance from the vertex at au scale and index 
a, which we also fit as a free parametre, is changed from 
-2 to 2. The best fit solution is looked for minimizing the 
\ 2 given by the sum of the square differences between the 
model and measured relati ve motions. We calculate t he \ 2 
from the equation (A. 2) in lGoddi fc MoscadellH i|2006h . 

The best fit solution is found to have a large open- 
ing angle of 54°. The cone axis is inclined at 37° from the 
line of sight and its position angle is 60° east from north. 
It is worth noting that this NW-SE axis coincides very 
well with the large scale jet seen in 2.1 /im hydrogen emis- 
sion IjAnandarao et ahl 120041 ). The position of the vertex is 
(a,<5)j2ooo.o = (06 h 09 m 06 s .9635, +21°50'41".3872). This po- 
sition has an offset of 150 mas from the peak position of the 
VLA continuum source to SW direction, but barely within 
the error of the absolute coordinates. This position is com- 
pletely inside the continuum emission. The best fit a is -1.7, 
but, this may be a mean value of each location, because R 
in equation (2) is also a function of the density of parental 
clump. Hence we need more detailed model, which includes 
flow-geometry and density gradient of clump, to determine 
the relation between R and R. 

This large opening angle solution and axis geometry can 
be interpreted by a dual structure of highly collimated jet 
and massive, large opening angle molecular outflow. Each 
types of flow is often ob served in massive sta r forming re- 
gion (e.g., the former in iBeuther et al.ll2002l ; the latter in 
IShepherd. Testi fc Starkll2003l ). and someti mes, they actu- 
ally c o -exist associated with single MYSO (jShepherd et al.l 
1 1999! ) . IBeuther fc Shepherd! i|2005h suggest that this co- 
existance can arise among the intermidiate phase on an 
evolutionary stage of early B-type star or B-type stage of 
mass accreting early O-type star. Recent MHD simulations 
also show that such a dual structure can be formed by a 
combin ation of acrretion and deflection of infalling m at- 
terfe.g jLerv et al.ll2002l ; ICombet. Lerv fc Murphy|2006l ) . In 
this case, the maser emission is more likely to be associated 
with the molecular outflow, since the jet is too collimated to 
excite the whole of the widely distributed maser features. In 
addition to this, the density of such a jet is significantly lower 
than that of an outflow in these MHD simulations, hence, 
massive outflow is more suitable than a jet for a maser ex- 
citation in terms of mass supply which supports to produce 
extremely dense gas layer widely. 

If we assume that observed deceleration of the V-wing 
is caused by swept-up of the cloud, density of outflow is 
estimated from equation (2) as 

"-How = ((-R) 2 + RR)n cloud Vfi 2 , 



1.4 x 10 10 (t^t)- 2 [cm" 3 ] 
kms -1 



(5) 



We gave n c i oud =10 7 cm -3 , R = 90 km s _1 , R = - 
1.5xl0 -6 kms -2 , R — 30 au, where n c i ou( j is a lower limit 
to the pre-shock density and others are given by the obser- 
vations in the case of IRAS 06061+2151. R is a half value 
of the largest relative proper motion in the V-wing and R is 
also observed deceleration of the V-wing. Compressed length 



R is a typical length of the path in which a maser feature 
moved during the observations. This length is slightly short 
to achieve a typical maser cloud density (10 9 cm -3 ) with 
the compression only, but if we take into account a mass 
conservation of the massive outflow, sufficient mass is pro- 
vided for a gas layer. Equation (5) indicates that if molecular 
flow has a typical velocity of a f ew tens of km s -1 (e.g ., 
IShepherd fc Churchwelll Il996al lbl; IChurchwell et all I2002T ). 
then its density is more than 10 6 cm -3 . This is the most 
dense case in the observed massive outflows. Such a molec- 
ular f low is, however , has not been detected in this region, 
yet (|Snell e t al. 1988]). This may mean the driving source of 
molecular flow is still compact and deeply embedded in a 
dense clump, and we can only see an collimated jet which 
can penetrates surrounding cloud. 

lAnandarao et all (|2004l ) regards their S4 source as the 
driving source of the jet, because it is thought to be a 
Class I B-type protostar which is thoght to form protostel- 
lar jet intrinsically and it is also located at the middle of 
the two knots. Our model fitting, however, suggests that 
the root of these jet/outflow system is located within the 
small UCHII region, G188.794+ 1 031 (see fig 1). I n fact , 
two .Ks-band sources, SI and S4 in lAnandarao et alj (|2004T ) . 
and G188. 794+1.03 1 are almost in alignmen t with the two 
knot s (see figur e 4 in Anandarao et all (|2004l ) and also figure 
4 in iBik et all (|2006T )). But if collimated jet and massive, 
embeded outflow actually co-exist, S4 seems to be inade- 
quate for their driving source. Because S4 can be seen in 
near- infrared (A"s-band) emission, it should not be deeply 
embeded. G188. 794+1. 031 is, on the other hand, only seen 
in 3.6 cm continuum and located at the blank of ifs-band 
emission. This indicates that G188. 794+1. 031 is still deeply 
embedded and not seen in near-infrared wave-length. This 
is consistent with our model mentioned above, and hence, 
the ionization source of G188. 794+1. 031 is more likely to be 
the driving source of the jet/outflow system. 

We estimate the dynamical time scale of the jet which 
is a rough measure of the lower limit on the age of driving 
source. The jet is traced by two, NW-SE, knots seen in 2.1 
/xm hydrog en emission which is a good tracer for moder- 
ate shocks (|Shull fc Beckwith|[l982l ). We assume a moderate 
shock velocity (10-30 kms -1 ) as a lower limit on the jet ve- 
locity, because this is only a terminal velocity after the de- 
celeration on t heir travel of 0.5 pc (from end to end) in dense 
molecular gas (jAnandarao et al.ll2004r ). We also assume the 
largest proper motion of maser features as an initial velocity 
of the jet, because it is seemed to reflect the molecular out- 
flow velocity whi ch is usually slower than the co -existing col- 
limated jet (e.g.. lLerv. Henriksen fc Fiegelll999h . With these 
assumption, estimatied time scale is ranged from 3.0xl0 3 
to 1.5 xlO 4 yr. The upper limit is not so small value, but it 
is still smaller than typical lifetime of ultra-compact phase 
of an ionized regi on around a massive star (~10 5 yr ; e.g., 
iHoare et ai1l2007h . 

4.3 Comparison of the objects on the jet axis 

The geometry in IRAS 06061+2151 suggests that the three 
sources on the jet axis (G188.794+1.031, SI and S4) may 
have some dynamical relation through the jet/outflow sys- 
tem. We cannot estimate the age of these three sources di- 
rectly, but an absence of associated ionized region is rough 
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measure of their growth. In addition to this, comparison of 
a spectral type between these sources give some decision on 
their formation sequence. The spectral type of the ionizasion 
source of G188. 794+1. 031 can be estimated from the excita- 
tion parameter which is calc ulated with 8.4 GH z continuum 
emission and equation (3) in lKurtz et al.l l| 19941 ). Calculated 
parameter is described as, 

where D is a source distance from us and we assumed the 
electron temperature as 10 4 K. An excitation parameter of 
1.92 correspo nds to that for a little more massive star than 
B3-type star (|Panagialll973r ). This suggests that the ioniza- 
tion source of this small UC H II region has lower limit mass 
that can produce a ionized region (7-8 Mq). 

In contrast to G188. 794+1. 031, other two Tfs-band 
sourc es have no overlayed radio counterpart (|Bik et al.l 
2006), and hence, no detectable UCHII region is associ- 
ated with them. Nevertheless, the spectral type of these two 
ifs-band so urces, which is estimated from colour-magnitude 
diagram by lAnandarao et al" | (|2004T ). is B1-B0 type which 
is enough early to produce an ionized region. This indicates 
that these two Ks-band sources are more massive than the 
ionization source of G188. 794+1. 031 which have lower limit 
mass for an ionization. If these spectral-type estimation is 
correct, the absence of any ionized region around these Bl- 
B0 type sources suggests that they are younger than the 
ionization source of G188. 794+1. 031, because they should 
evolve more rapidly than, or at least as rapid as, B3-type 
source. Taking into account this result, the ionization source 
of G188. 794+1. 031 seems to be born at first in these three 
sources. In this case, if the ionizing source actually drives 
the jet, the formation of two B-type protostars, especially 
S4 source which adjoins G188. 794+1. 031, may be affected by 
the jet /out flow system. In this case, it can be a good example 
of sequential star formation in a massive cluster. In order to 
prove such a sequential process in IRAS 06061+2151 quan- 
titatively, we have to determine their geometry more accu- 
rately, in particular, the shape and the size of the UCHII 
region and the absolute location of the maser. 



5 CONCLUSION 

We made VLBI observations of H2O masers in massive 
star forming region IRAS 06061+2151 with JVN from 2005 
May to 2007 October. The H2O masers are clustered within 
l"xl" area around extremely young UCHII region. They 
are divided into several groups based on their radial veloci- 
ties, and show expanding motions. Their alignment and de- 
tected three-dimensional motion indicate they are associated 
with a bipolar outflow from the source in the UC H II region. 
The southern blue-shifted group shows an arc-like alignment 
within several tens of au. Their continuous alignment indi- 
cate they are excited in a single gas layer. They also show 
two clearly different kinematics, which are the vertical part 
(V-wing) of 90-180 kms -1 and the horizontal part (H-wing) 
of 30-65 kms" 1 . This may be caused by far small (several 
tens of au) and clumpy (10 7 -10 8 cm -3 ) structure of parental 
cloud core. Such a clumpy structure is expected to be in- 
duced by strong instability of shock front (e.g., expansion of 



UC H II region) . Further high-resolution MHD simulation is 
desirable to examine this effect. 

The biconical flow has an opening angle of 54° and 
the axis along NW-SE direction based on our model fitting. 
This NW-SE axis coincides very well with the large scale jet 
seen in 2.1 /^m hydrogen emission. These geometry and large 
opening angle are thought to indicate the dual structure of 
highly collimated jets and massive and largely opened molec- 
ular flow in this region. We estimate density of the outflow 
from the deceleration of the V-wing of ~10 6 cm -3 , although 
such a massive outflow has not been detected in this region 
yet. This may mean such a flow is still compact and deeply 
embedded in the dense clump, and we can only see an col- 
limated jet which can penetrates surrounding cloud. In this 
point of view, we conclude that these two types of flow is as- 
sociated with the ionizing source of G188. 794+1. 031. There 
are, in addition, two B-type protostars located on the jet 
axis. These two sources are seemed to be younger than the 
ionizing source from the absence of any ionized region and 
comparison of a spectral type. This time sequance and their 
geometry may indicate that the formation of two protostars 
is affected by the jet and outflow. With additional obser- 
vations which has higher resolution and accuracy, they can 
be a good example of sequential star formation in a massive 
cluster. 
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